WELL-POSEDNESS AND AVERAGING OF NLS WITH 
TIME-PERIODIC DISPERSION MANAGEMENT 



PAOLO ANTONELLI, JEAN-CLAUDE SAUT, AND CHRISTOF SPARBER 



Abstract. We consider the Cauchy problem for dispersion managed nonlinear 
Schrodinger equations, where the dispersion map is assumed to be periodic 
and piecewise constant in time. We establish local and global well-posedness 
results and the possibility of finite time blow-up. In addition, we shall study 
the scaling limit of fast dispersion management and establish convergence to 
an effective model with averaged dispersion. 



(1.2) 7W 



1. Introduction 

In this work, we study the Cauchy problem for the following class of dispersion 
managed nonlinear Schrodinger equations (NLS): 

(1.1) idtu + -f{t)Au+\u\P-\^0 u{tQ,x) = ip{x). 

for {t, x) G X R'', a given to G IR.+ , and p > 1 to be specified below. In addition, 
we assume that the dispersion map 7 is 1— periodic, i.e. j{t 4- 1) = 7(0' 
piecewise constant: 

'7+ 0<ti^t+, 
-7- t+<ti^l, 

where 7± > are some positive constants and t+ G (0, 1). 

Our main motivation stems from models in nonlinear fiber optics. Indeed, the 
case of a cubic NLS, i.e. p = 3, in d — 1 spatial dimension naturally arises as an en- 
velope equation for electromagnetic wave propagation in optical fibers exhibiting a 
(weak) Kerr nonlinearity. In this context, the variable t G M-i- actually corresponds 
to the distance along the fiber, and a: G K. denotes the (retarded) time. The coeffi- 
cient 7(<) consequently models a periodically varying dispersion along the fiber, cf. 
[2, 32]. Dispersion managed NLS in d = 2 spatial dimensions are also physically 
relevant, see, e.g., [1]. The technique of dispersion management was invented to 
balance the effects of nonlinearity and dispersion in such a way that stable nonlin- 
ear pulses (solitary waves) are supported over long distances, cf. [22, 25, 26, 34]. 
Due to the enormous practical implications, there is a huge literature concerned 
with the qualitative properties of (1.1). Most often, however, the results arc based 
on non-rigorous asymptotics and/or numerical simulations. A notable exception is 
the regime of so-called strong dispersion management where several rigorous results 
are available for the corresponding asymptotic model, see the discussion below. 

In contrast to that, we shall work directly on the dispersion managed NLS and 
in the following prove several rigorous results concerning the well-posedness of (1.1) 
and its asymptotic behavior in the case of rapidly varying dispersion. Having in 
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mind the physics background of fiber optics, we shaU focus on the well-posedness 
theory for (large) data </? € L^(R'') or </? € H'^{R'^). By multiplying (1.1) with u, 
integrating w.r.t. x £ M.'^, and taking the imaginary part of the resulting identity, 
we (formally) obtain mass conservation: 

•)llL2(ffid) = ||</3||i2(Rci). 

Note that in the case without dispersion management, i.e. j{t) = 7 G M, equation 
(1.1) also conserves the energy 

(1.3) E{t)^- I \Vu{t,x)\'^dx 3_ / \u{t,x)Y'+^dx^ E{Q). 

Thus, in the defocusing case 7 < and when p < ^35 , one immediately infers a 
uniform bound on the norm of u, prohibiting the appearance of finite-time blow- 
up. In view of (1.2), we can think of (1.1) as switching in between focusing and 
defocusing behavior (with respective dispersion coefficient 7+ > and — 7_ < 0) 
and we consequently expect the appearance of finite time blow-up. That this is 
indeed the case will be proved in Section 3.2. 

Remark 1.1. Note that this switching between focusing and defocusing behavior 
makes our problem very different from the NLS type models studied in [16, 23]. 
In [16] the author considers an equation of the form (1.1) with time-dependent 
dispersion j(t) ^ 0, possibly vanishing (with finite order) on a discrete set of 
points, a typical example being j{t) = |t — ^oj^ for some A > 0. In [23] an NLS type 
model with time-dependent coefficients but with an additional constant coefficient 
third order spatial derivative is studied. The behavior of this model is similar to 
the KdV equation. 

Even if finite time blow-up in general can not be ruled out, one may still wonder 
if sufficiently fast switching between the focus and defocusing step can at least 
delay, the appearance of blow-up. In order to gain more insight we shall therefore 
study the scaling limit corresponding to fast dispersion management, i.e. we shall 
consider 

(1.4) tdtu,+-f(^^^Au, + \u,\P-^Ue = 0, u,{to,x) = ip{x). 

where < e ^ 1 denotes a small parameter. This regime has been studied using 
formal asymptotics and numerical simulations in, e.g., [6, 7, 35]. As £ — > 0+ we 
expect the behavior of the solution u'^ to be close (in some sense to be made precise, 
see Section 4) to the solution of the averaged NLS 

(1.5) idtuo {j)Auo + \uo\^~^uo 0, uo{to,x) ip{x), 
where we denote by 

(1.6) (7) := / 7(T)dT, 

Jo 

the average dispersion coefficient (which can be either positive or negative). In 
the case of mean zero dispersion (7) = this scaling limit provides a possible 
explanation for the stabilizing effects of dispersion management, see Corollary 4.6. 

The situation above should, however, be distinguished from the case of strong 
dispersion management, cf. [25, 17] for some physical motivation. One thereby 
considers a dispersion map of the form j-7(|), which, as e — > 0+, leads to an 
effective description by a non-local equation, originally introduced in [17]. This 
model has been rigorously studied by several authors, see, e.g., [14, 20, 30, 31, 34]. In 
particular, it provides a mathematical basis for the definition of dispersion managed 
solitons [14]. 
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Remark 1.2. A similar situation is analyzed in [5, 13], where the authors con- 
sider NLS with fast random dispersion management. More precisely, they consider 
dispersion maps of the form ^-7 (^), where 7 is is a (smooth) stationary random 
process, and prove convergence of the equation towards an NLS with white noise 
dispersion. Finally, we want to mention that in [11, 12] the somewhat dual prob- 
lem of NLS with (rapidly) time-oscillating nonlinearity has been analyzed (sec also 
[9, 28] for a related studies on the KdV equation). 

This paper is organized as follows: In Section 2 we shall collect some basic 
properties of the linear equation to be used in the following. Section 3 establishes 
global well-posedness in the sub-critical regime. In addition, we set up a local 
well-posedness result in theory and prove the existence of finite time blow- 
up. Section 4 is concerned with the fast dispersion limit. Finally, we shall collect 
some concluding remarks on possible generalizations and closely related problems 
in Section 5. 

2. Basic properties of the linear equation 

2.1. The linear propagator. Before studying the nonlinear Cauchy problem we 
shall collect some basic facts about the associated linear equation 

(2.1) idtUlin + 7(t)Auiin = Ulin(to, x) = ip{x). 

In the following, we shall denote by 

r(t,s) / 7(T)dT 



the cumulative dispersion on the timc-intcrval [s,t] C M and the associated propa- 
gator U{t,s) by 

(2.2) C/(i,s)/(x) :=e^i^(*'^)^/(x)= / e-'i^(*'^)l«l7(0e"-«dC, 



where / denotes the Fourier transform of / G L^(R'^). We directly infer the follow- 
ing result. 

Lemma 2.1. Let to (zR be fixed. Then the mapping t i— > U{t,to) defines a family 
of strongly continuous unitary operators on L^(R''), such that for all t G M_|- it holds 

Ulin{t,x) = U{t,to)ip{x). 

Explicitly, we find 

provided T{t,tf)) 7^ 0. 

Some remarks are in order: First, one should note that for all s, t G M: 

U{tM)U{tQ,s) = U{t,s), 

and we also have 

U{t,s)* = U{t,s)-^ = U{s,t). 
However, U{t, s) 7^ U{t — s, 0), since 

l^t pt — s nt — S 

j{T)dT = / 7(t + s)dT 7^ / j{T)dT, 

Jq J a 

unless 7(t) — const, for all t G M. In other words, U{t,s) is not a group. Second, 
we note that for ip G L^{R'^), it holds 

dtue BV{R;H-'^{R'^)), 
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in view of (1.1) and the fact that "f{t) is piecewise constant (including jump dis- 
continuities). In particular, dtu is not continuous in time, in contrast to the case 
with constant dispersion. 

Next, we recall the decomposition 

7(i) = (7)+7o(t), 

where (7) G R denotes the average defined in (1.6), and 7o(0 is 1— periodic with 
mean zero. Using this, we can write 

(2.4) r(t,s) = (7)(t-.s)+ / 7o(T)dr. 



We consequently infer that if (7) = 0, then r{t,to) = for all t € M such that 
(t — to) G N. On the other hand, if (7) ^ 0, then r{t,to) vanishes at most finitely 
many times, since 

(2.5) -(7- + (7»«; / 7o(T)dT s;7+- (7), 

J to 

and thus |r(t,<o)| ^ |^|, provided t ^ to is sufficiently large. Also note that if 
to € N, i.e. if we start our time-evolution at the beginning of focusing step, then 

d{t) := f 7o(r)dr ^ 0, 

Jto 

for all times t G M, whereas if we start at t = to at the beginning of a defocusing 
step, then ^{t) ^ for all times. 

Remark 2.2. In the particular situation of mean-zero dispersion, i.e. (7) = 0, the 
solution liiin is found to be 1— periodic in time. Indeed, in view of (2.2), we infer 
that the Fourier transformed solution is given by 

t 

2 



uiin(t,0 = exp 7o(r)drj ^(0, 

which satisfies uiin{n,x) ~ '^{x), for all ri G N. 

2.2. Dispersive properties. From what is said above it is clear that [/(t, s) in 
general does not allow to infer uniform (in time) dispersive estimates (analogous to 
the usual Schrodinger group), since for arbitrary s,i G R 

t 

' t- s\ 



\T{t,s)\= J ^{T)dT 

due to fact that ^{t) changes sign. 

Remark 2.3. Clearly, this would not be the case if \^{t)\ > and 1— periodic, in 
which case the behavior of (1.1), for all times t G K, would be either focusing or 
defocusing (but with time-dependent dispersion coefficient), see, e.g., [16]. 

The fact that 7(t) is assumed to be piecewise constant allows us to obtain the 
following result. 

Lemma 2.4. Let t^ G (0,1) be given and t, s d {n,n + 1^], or t,s ^ (71 + 71 + 1], 
for n G N. Then, for t ^ s, it holds 

\\U{t,s)f\\L^(^.)<\t-s\-^''~\\f\\L.(^.y 

The point is that both t and s have to be in the same time-interval corresponding 
to either a focusing or defocusing step. 
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Proof. The proof follows directly from the representation formula (2.3) and the fact 
that if t, s G (n, n + t+J, or t, s G (n + t-^,n + 1] we have 

r(i,s) = ^ 7(T)dT = ±7±(t-s), 

in view of (1.2). □ 

Clearly, we have that 
(2.6) [j {in,n + t+]U{n + t+,n + l]) =R, 

riGZ 

and thus we can split the time-axis R into a union of sub-intervals on which U(t, s) 
allows for the usual dispersive behavior of the Schrodinger group. As a consequence, 
we know that the usual Strichartz estimates hold on each such time-interval [19, 21]. 
To this end; let us recall that a pair {q, r) is admissible if 2 ^ r ^ (resp. 
2 < r oo if d = 1, 2 ^ r < oo if d = 2) and [10]: 



q \2 r 

Lemma 2.5. Let rt G N and t+ G (0, 1) be given. Then, for each admissible {q,r), 
((7i,ri), {q2,T2), and each time-interval /„ C {n,n + i+], or /„ C (?i + n + 1], 
we have: 

1. There exists Cr — C(r, /„), such that for to G /„ and for any f G L^(R'*) it holds 

2. There exists Cri,r2 = C{ri,r2,I), such that for any F G L"^^ (I„; L"^^) it holds 
[ U{t,s)F{s)ds < C'.^.JI^^II 



L«2(7„;L'^2 (Erf)) 



Again, we can only conclude the existence of Strichartz estimate locally on each 
time-interval corresponding to, either, a focusing or defocusing step. 

Remark 2.6. One may want to compare this to the case with .smooth periodic 
dispersion management. For example, let jit) = cos(i) and to = 0. In this case we 
find 



|4^sin(t)|'i/2 7R. 

for which the corresponding dispersive estimate is again valid only for small times. 



I.e. 



||;7(t,0)/||i^(K.) < |t|-''/'!|/||Li(Erf), for |t| < |. 

This situation is very similar to the case of a Schrodinger equation with quadratic 
confinement [8] . The corresponding Schrodinger group is formally given by 

Sit, 0)ip{x) = e''*(^+l"l'V(^), t e R+, 

the kernel of which can be explicitly computed via Mehler's formula. Due to the 
existence of eigenfunctions, this time-evolution enjoys dispersive properties only for 
sufficiently small < |t| < S, sec [8] for more details. 
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3. Well-posedness results 

3.1. The sub-critical case. As we have seen, the propagator U{t,s) does not 
allow for uniform dispersive estimates. Nevertheless, one can prove global well- 
posedness of (1.1) in the subcritical case. 

Theorem 3.1. Consider (1.1) with ip e L^(M'*) and 1 < p < 1 + |. Then, for any 
to G K., there exists a unique solution u e C(R, L^(R'^)) satisfying 

\\u{t, ■)llL2(M'i) = ll</'llL2(Rd). 

In addition, for each admissible pair {q, r) and for any compact time-interval I C M.' 

!l"llL|Lj(/xR<i) ||uo||i2(]B.d)). 

Proof. Using Duhamel's formula the solution of (1.1) can be written as 
(3.1) u{t,x) ^U{t,to)ifi{x) +i I U{t,s)\u{s,x)Y'-^u{s,x)ds. 



Next, we use the decomposition (2.6) to split the time-interval [tQ,t] e M into 
countably many sub-intervals 



[io,i]- U ([to,i]n/i)u([to,t]n/2)), 



with l}^ = {n,n + t+] and if^ = {n + t+,n + 1]. In each of the intervals 

we are able to apply the Strichartz estimates stated in Lemma 2.5. Indeed, for 

1 < p < 1 + I , the pair 

('^'")^(d(^''^. 
is admissible and hence we have 

. 4p ^ I1"o||l2 + |||u|^"^"||LiL2(/xRci) 

< ||^/,o|U. + |/r||^/,r 4. 

where a = i+d-dp ^ q jjj,^^^,^ if / c I^'^ is sufficiently small, a standard continuity 
argument implies that there exists a solution u to (1.1) in / x such that 

||u|| 4p s$ C(||mo||l2(r<'))- 

The right hand side only depends on the L^(IR'')-norm of the solution, which is 
uniformly bounded for all times (and in fact equal to ||(y9||/^2). Continuity therefore 
implies that on each subinterval the solution exists for all times t in one of the 
sub- interval and, in addition, we have 

By considering the union of sub-intervals /„ we consequently infer the existence of 
a solution u e L°°(]R; L2(M")). 

In order to obtain the asserted continuity in time, we consider two different times 
t2 7^ ti G M.^ for which we need to show that 

||u(ti, •) - "(^2, •)l|L2(Rn) ^ 0, as |ti-t2|->0. 

In view of Duhamel's formula (3.1) this requires 

\\Uit,,to)^ - U{t2,toML2(M^) 0, 
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which is nothing but the strong continuity of U{t, to) stated in Lemma 2.1. For the 
second term on the right hand side of (3.1) we write 

ti rt2 

U{ti, s)\u{s, x)\^~^u{s, x) ds — / U{t2,s)\u{s,x)\^~^u{s,x)ds — ii + 12, 

3 J to 

where 



ii = (?7(ii,io) - C/(t2,io)) / U{to,s)\u{s,x)\P-^u{s,x)ds 



to 



i2= U{t,s)\u{s,x)\P-^u{s,x)ds. 



and 



By using the strong continuity of U (t, s) for ii, and Strichartz estimates for i2 (after 
sphtting the time-interval [^2,^1] into unions of and /^), we see that both are o(l) 
as 1^1 —t2\ 0. Thus, the solution is continuous in time with values in L^(M''). □ 

3.2. Local well-posedness in and finite time blow-up. In order to allow 
for p = 1 + 2, which, in particular, includes the physically relevant situation of 
a cubic nonlinearity in d = 2 spatial dimensions, we require tp € i7^(M'*). In the 
following, we shall always denote by (a)+ the positive part of a e M. 

Lemma 3.2. Consider (1.1) with (p e H^{M.'^) and 1 < p < ^ + jj-^j^ ■ Then there 

exists a T* = T*{\\(f\\H'^) > to and a unique solution u £ C{[to,T],H^(R'^)), for all 
to < T < T* . The solution is maximal in the sense that if T* < +00 then 

^hm \\Vu{t, •)||L2(Rd) = +00. 
Note that in dimensions = 1, 2 we can allow for any exponent 1 < p < -|-c». 

Proof. Differentiating (1.1) with respect to a; we sec that Vu solves the following 
integral equation 

Vu{t,x) = U{t,to)\/ip{x) +i f U{t,s)G{u{s,x))ds, 

J to 

where 

G{u) = ?-^\u\P-^Vu + ?-^\u\P-^u'^Vu. 

The latter satisfies 

G(w) < \u\P-^\\/u\. 

By decomposing, as before [<o,^] hito unions of and and applying Strichartz 
estimates together with Soblev's imbedding on each of this sub-intervals, yields the 
existence of a local in-time solution u{t, •) S H^{M.'^). The existence T* > time 
thereby depends on □ 

Clearly, if at the initial time to G K it holds: 7(^0) = ^ 7-, i-c if we start with 
a defocusing step, then we know that the solution exists at least up to T* ^ 
This follows from the fact that j{t) is constant (and negative) for all t £ [to,t+] 
and thus the energy (1.3) is conserved on this time interval which in turn provides 
a uniform bound on the norm of u. 

Remark 3.3. In the case of constant dispersion and if p < 1 + ^, one can infer 
T* = -|-oo, i.e. a global solution, using Strichartz estimates and the conservation 
of energy. Since conservation of energy does not hold in our case one has to argue 
differently. Namely, ii (p e i?^(M'') and Hp < 1 + |, the global solution 
constructed in Theorem 3.1 admits propagation of regularity and thus we indeed 
obtain that u e C(E, i?i(R'')). The proof of this result is given in [10, Theorem 
5.2.11. 
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In the following we denote by Q G ^^(M'^) the unique radiahy symmetric solution 
of the following elliptic equation: 

(3.2) AQ-O + Qi+^=0. 

The existence of such ground states Q has been proved in [3, 4]. 
Theorem 3.4. Letp=l + ^ and ip e H^{W^) with 

Ml^<i+'"'\\Q\\l-- 

Then T* = +oo. Moreover, this criterion is sharp. 

Up to the inclusion of the scaling factor 7^ , this criterion is analogous to 
the case without dispersion management, cf. [24]. We see that the stronger the 
dispersion in the focusing step, the larger the mass of the initial data can be in 
order to infer global existence. 

Proof. With the local in time existence result for solution at hand, we need 
to show that |j Vu(t, •)||i2(-][jd') < +00 for any finite time t S R. To this end, let 
us assume for the moment that at t = to we start with a defocusing step, i.e. 
7(^0) = ^7- < 0. Denoting by [iol G N the ceiling of to G K, we consequently have 
that for all t e [to, \to]] the following energy is conserved: 

E^{t):^\[ \Wu{t,x)\^dx + -^ [ \u{t,x)\'+Ux = E^{to), 

where we have set p = 1 + A/d. Next, let us recall the Gagliardo-Nirenbcrg type 
inequality (see [33]): 

(3.3) WfWlttit ^ (1 + IIOilZ'^'ll/!l'/^'ilv/lli., 

where Q is the solution of (3.2). Using that u(to) G i7^(]R'*), by assumption, 
together with the fact that for all times < G M: ||u(t)||i2 = j|</3j|L2 < +00, the 
inequality (3.3) implies that E-{t) < +00 for all t € [to, [to]]. We consequently 
obtain that u(t) e H^iW^) exists up to t = [to] € N and that 

(3.4) \\\/u{\to^,■)\\h<E^ito), \\ui\to^,■)\\ltti%<E-{to)■ 

On the time interval [[to], [to] +t+] we are in the focusing regime. The associated 
energy is 

E+it)-^^ [ |Vw(t,a;)|2dx-— ^ / |u(t,x)|2+^^dx = £;+([to]) < +00, 

where the last inequality follows from (3.4). Using again (3.3) we infer that for all 
t e [[to], [to] +t+] it holds: 

i?+([to]) = E+{t) l|V^.(t, .)||i2. 

This consequently implies that if < 7+'^^^||Q||l2 we have 

/ II ||4/d\-l 

||V^.(t)||i2^2 7+-^^ i?+([to])<+oo, 

for all t G [[to] , [to] + 1_(.]. Thus, we infer the existence of u{t) E H^{M.'^) up to the 
time t = [to] + t+ after which we are again in the defocusing regime and the same 
argument as above applies. Clearly, this also shows that we could have started at 
time t = to with a focusing step, instead of a defocusing. In both cases we obtain 
that ||VM(t)||^2(-jj£i) < +00 for any finite time t € R. 
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Next, let us show that this result is indeed sharp, i.e. we need to find an initial 
data Lp S H^{W^) such that ||iy3||^2 = 7+''^'^||<3||l2 and the corresponding maximal 
time of existence for u is T* < +cx). By the pseudo-conformal symmetry of the 
mass-critical NLS with p = 1 + | it is well-known that (see, e.g., [10, 24]): 

(3.5) Va{t,x)^{l-atyiQ 
is a solution to the mass critical focusing NLS 

4 

idtVa -f 7+AWa +\Va\'^Va= 0, 

with initial datum 

satisfying ||wa(0, OIU^ = I+^^'^WQWl^- In addition, we see from (3.5) that Va blows 
up at time t = -. Now if to is such that 7 (to) = 7-i-j i-C- we start with a focusing 
step, then all we need to do is choose a = a* such that 

1 r 1 

— <t+- {to}, 
a* 

where [0, 1) 3 {to} = to~ L^oJ is the fractional part of to. Setting u(to, x) = Va, (0, x) 
ensures that the solution to (1.1) will blow up before the time [to J +t+, i.e. before 
switching to the defocusing regime. If on the other hand we start in a defocusing 
step 7(to) = —7-, then we choose u{tQ,x) = u)(to,x), where w solves the time- 
reversed defocusing NLS 

4 

idtw + 7_ Au; — = 0, 

with initial data w([to] , x) = Va, (0, x). □ 

In particular, the proof shows that a defocusing step in general can not prevent 
the appearance of finite time blow-up. Note, however, that the construction of 
blow-up solutions given above is based on the pseudo-conformal symmetry of the 
mass critical NLS and the associated blow-up scenario is known to be unstable, see, 
e.g., [24] for further discussions. 

Remark 3.5. Also note that in the proof Theorem 3.4 we have used the fact that 
7(t) is piecewise constant several times. Indeed, the situation where 7(t) is smooth 
and 1— periodic would be considerably more complicated. 

4. Averaging for fast dispersion management 

4.1. Preliminiaries. In this section we consider the £— scaled NLS (1.4). Clearly, 
we have that for every fixed e > all the results of the foregoing section remain 
valid. In the following, though, we shall consider the scaling limit as e — >■ 0+, i.e. 
the regime of rapidly varying dispersion. 

As a first step, we shall show that the expected averaging result is true for the 
linear equation. To this end, we denote by 

rt/e 

[/,(t,.s) = e*r=(*.^)A, r.(t,s)=e/ 7(T)dT, 

Js/e 

the propagator associated to the fast-dispersion map 7 (7). In addition, 

C/o(t, s) = e^^-'X^-^)^ ^ e"<^>^ o e-*^<^>^, 

is the propagator associated to linear Schrodinger equation with averaged dispersion 
(7) e K defined in (1.6). Note that in fact Uo{t, s) = Uo{t — s, 0) and thus a group. 
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Lemma 4.1. For any e M and f G H'^{R'') with a e N we have 

lim sup \\Ue{t, s)f - Uo{t, s)/||H-(Md) = 0. 
"^^0+ t,seR 

Proof. Using (2.4) we can decompose as 

ft/e 

T,{t,s) = {j)it-s) + e / 7o(r)dr = (7)(t-s)+ei?,(<,s), 

where we note that s) G L°°(M'') uniformly, in view of (2.5). Fourier transfor- 
mation and Plancherel's identity allow us to write 

2 _ 



sup \mt,s)f-Uo{t,s)f\\j^ 

t,s£S. 



2 



= sup / (l + l^l^)- e'Wtt-«JKI g^e#e(t,s)_l \fiO\ d^- 

The assertion then follows by the Lebesgue dominated convergence theorem and 
the fact that 

lim (e'''^^^*'''^ -l) =0, 

E-i-0+ V / 

pointwise for all t,s E M.. □ 

To prove the desired averaging result as e — >■ 0+, we will require sufficiently 
smooth solutions Ue{t,-) G H'^{W^) with cr > |. In the following we shall concen- 
trate on the physically relevant cases of d ^ 3, spatial dimensions. The generaliza- 
tion to higher order dimensions will be indicated below. For d ^ 3 it is sufficient 
to consider solutions in iJ^(IR'^), whose existence is guaranteed by the following 
lemma. 

Lemma 4.2. Let (p e iJ^(M^) and 2 ^ p < oo. Then, for any e > Q, there exists 
a time T* = T*(||(/3||//2; e) > tg and a unique solution Ug € C{[to,T], H'^{M.'^)), for 
all to < T < T* . The solution is maximal in the sense that if T* < +oo then 



Furthermore 



lim ||ue(i, 



inf T*(||^||ff.;e)>0. 



Proof. The proof follows from arguments analogous to those given in the proof of 
Theorem 1.4 in Chapter 6 of [29] (see also [10, Section 4.8]). In particular, this 
establishes the fact that t > 0, bearing in mind that supp^^^j^ \\Ue{t, s)\\ = 1. □ 

Standard arguments also imply that for (p £ H'^(R.'^) there exists a maximal time 
of existence T, = r*(||(^|| j/2) > to and a unique solution uo € C{[to,T], H^(R'^)), 
for all to < T < T^., satisfying the averaged equation (1.5), cf. [10, 24] for more 
details. Depending on the power of the nonlinearity p and on the sign of (7) G M. 
the maximal time of existence for the averaged equation might be infinite or not. 
In particular, for ff-'^-subcritical nonlinearities in the defocusing case (7) < we 
have = +00, cf. the results given in Section 5.3 of [10]. 

4.2. Averaging of NLS with rapidly varying dispersion. The main result of 
this section is as follows: 

Theorem 4.3. Let d ^ 3, ip e H'^{R'^), and 2 p < cxj. Denote by uo S 
C([to, r*), -ff^(M'^)) the maximal solution to the averaged equation (1.5). Ln addi- 
tion, let e C{[to,T*),H'^{W^y) he the maximal solution of (1.4) for given e > 0. 
Then, we have that 

u/^uo, mL'=^{\to,T],H\R'')), 
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for all to <T <T^. 

In particular, we know that as e — >■ 0+ the solution of the original equation 
(1.1) can not blow up before > 0, the maximal time of existence of a smooth 
solutions to the averaged equation. 

Proof. In the following, we denote the nonlinearity by f{u) = |it|P~^M, for simplicity. 
In view of Duhamcl's formula (3.1) we need to estimate the difference 

Ue{t, x) - Uo(i, x) = ([/e(t, to) " Uo{t, to))ip{x) 

+ i / Usit,s)f{us{s,x)) - Uoit,s)f{uois,x))ds. 

J to 

For the first term on the right hand side we can directly use Lemma 4.1, whereas 
the second term can be rewritten as 

/ Ue{t,s)f{us{s,x)) - UQ{t,s)f{uo{s,x)ds =If 

Jto 

where 

11= f U,{t,s){f{u,{s,x))-fiuo{s,x))) ds, 

Jto 

and 

I|= f iU,it,s)~Uoit,s))fiuo{s,x))ds. 

Jto 

We first start by deriving an estimate in L^(R''), i.e. a — 0. Using Minkowski's 
inequality, 

mhrLl < \T-to\\\mt,s)-Uo{t,s))f{uo{s))h^^Ll- 
As long as uo{t, •) e H^{«.'^) ^ L°°{M.'^) for d ^ 3 we have 

I1/('"o)||l2 ^ ||uo|Il^^||wo|Il2 < +00, 

and thus Lemma 4.1 implies that for e = s{T) sufficiently small there exists a 
(5(e) > such that 

Next, wc note that for all u, v e C it holds 

\f{u)^f{v)\^c{\ur^ + \vr')\u-v\. 

With this we can estimate Xf , using again Minkowski's inequality and the fact that 
Ue is unitary on i^, via 

\\m\Lf-Ll < ll/("e) - ./(mo)|1l;l2 < i\\Ue\\l}l + \\uo\\l^Jj\\u, - Uo\\l1l2 

As long as Ue,uo e H'^{R'^) ^ L°°(R^) this implies 

(4.1) \\u, - uo\\l^l2 < 5{e) + C\\u, - uoIIlil^- 

Next, we recall the following result proved in [11] (see also [9, 28]): 

Lemma 4.4 ([11, Lemma A.l]). Let T > and 1 ^ p < q ^ oo and A,B ^ 
some constants. Assume that f G L'^{0,T) satisfies, for all t G (0, T); 

ll/llL<-(0,t) ^A + B\\f\\L^i^o,t)- 

Then there exists a K = K{B,p,q,T) such that ||/||L<i(o,t) ^ AK . 
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Using this Lemma with q = oo and p = I, we infer from (4.1) that 

- uqWl^li < S{e) 0. 

To obtain the analogous estimate for H'^ with a = 1, 2 we differentiate the equation 
w.r.t. X G Mf^, using the fact that the nonhnearity is sufficiently smooth. Then, 
similar arguments as before imply that 

(4.2) -uo|U-H2 0. 

The result then follows from a continuity argument similar to the one given in 
[11]: Fix < r < and set M = supo^t^^- ||?io||H2. It follows from Lemma 4.2 
that for ||<p||//2 < M there exists a r > to such that for any e > 0, exists on 
[to , r] and 

(4-3) sup ll-UelUoc^t^^^j^^s) < ||(^||^2. 

0<e<l 

Next, let to < ^ ^ 2^ be such that exists on [to,^] for e sufficiently small and 

limsup\\ue\\L°-{{to,e),H2} < +oo- 

Note that r = £ is always a possible choice. Then we deduce from (4.2) above that 

from which we consequently infer that 

\\uS,-)\\h- ^M, 

for £ <C 1 sufficiently small. Applying Lemma 4.2 to the NLS (1.4) translated by £ 
we deduce that for e sufficiently small, exists on [to , r + £] and that 

sup ||We||ic»((t^^^+£).^2) < \\(p\\h^-. 
0<e^l 

In other words, the estimate (4.3) holds with £ replaced by £ + t, provided that 
£ + T ^ T < . Iterating this argument we infer that 

sup ||Ue||ioo((tj, ^.)^^2) < \\(p\\h^, 
0<e^l 

and the assertion is proved. □ 

As a first corollary we can state the following result valid in arbitrary spatial 
dimensions c? G N, provided the nonhnearity is sufficiently smooth. 

Corollary 4.5. Let d e N, f e H'^{R'^) with N B a > ^ and assume that the 
nonhnearity satisfies p ~ I + 2a with a € N. Then 

u/^^uo, mL°-{[to,T],H^m^)), 

for all to < T < T*, the maximal time of existence of solutions uo{t, •) G H'^{R'^), 
satisfying averaged equation. 

Proof. In higher dimensions, existence of smooth solutions uo,u^{t,-) G H'^{M.'^) 
with c > I can be proved as in [10, Section 4.10], provided p = 1 + 2q; with a G N. 
The averaging result then follows by the same arguments as given in the proof of 
Theorem 4.3 above. □ 

More interestingly, in the case of zero average dispersion (7) — the averaged 
equation (1.4) simplifies to an ordinary differential equation 

(4.4) idtuo + |uo|^"^uo = 0, uo(to, x) = '^{x), 
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which has been formally derived in [7] (including higher order corrections in e). 
Equation (4.4) can be solved explicitly, resulting in the following corollary (stated 
for d ^ 3, for simplicity). 

Corollary 4.6. Let d ^ 3 and (7) = 0. Under the same assumptions as in Theorem 
4-3 we have, for any compact time-interval I C K, that 

where 

wo(i,cc) = (^(x)e*(*~*°)l^(")l"^\ 

Proof. In the case (7) = 0, Theorem 4.3 implies that converges, for e — > 0+ to 
the solution of (4.4). Multiplying the latter by uq and taking the real part, we find 

i9t|uoP = 0, uo(io, x) = ip{x), 

and thus \uQ{t,x)\ = \'f{x)\ for all i S R. Writing uq — |tio|e*'' and, analogously, 
<p = |<p|e*^°, we find the following equation for the phase 

dte = \^{x)Y'-^e, 0{to,x)^0o{x). 

Integration w.r.t. t then yields the result. □ 

The particular form of uq found in the case of zero average dispersion corresponds 
to a solution which does not disperse. Indeed, the spatial density (corresponding 
to the energy density of an electromagnetic pulse) is seen to be time independent, 
i.e. |uo(t, a;)P = |(/3(a;)p, for all t G M. However, the solution uq oscillates with an 
x-dependent frequency u!{x) = |(y9(a;)|P~^. From the physics point of view. Corollary 
4.6 provides a possible justification for the stabilizing effect of dispersion manage- 
ment in optical fibers with mean zero dispersion. Note, however, that this effect 
should be distinguished from the ones established in [34] . 

Remark 4.7. The case for rapidly varying mean zero dispersion (7) ~ is radi- 
cally different from the corresponding situation found for NLS with time-periodic 
nonlinearity management [11, 12]. In the latter case, the effective model obtained 
after averaging is given by a linear free Schrodinger equation, whose solution is 
purely dispersive, in contrast to dispersion managed NLS. 

5. Concluding remarks 

5.1. The case \j(t)\ > 0. In this case, the behavior of (1.1) is purely focusing or 
dcfocusing. In particular, wc have that, for any given to G M, the mapping 

t ^ T(t,to) e M 

is strictly monotone and we can define a new unknown 

v(t,x) = u{T(t,to),x), 

which solves 

(5.1) idtv + Av + K{t)\v\P-^v = 0, v{0,x) = ip{x), 

with coefficient K{t) = Equation (5.1) is a NLS with time-dependent non- 

linearity management, similar to the models studied in [11, 16]. In particular, if 
7(i) < for all t G M, this equation is dcfocusing and one can prove global in time 
existence of solutions in H^{M.'^) along the lines of [16]. 

Of course, the (physically and mathematically) most interesting case of nonlin- 
earity management is the one with sign changing time-dependent coefficient K{t). 
Such a situation, however, is no longer equivalent to the one with dispersion man- 
agement. 
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5.2. Possible generalizations. Let us mention that it is straightforward to gen- 
erahze all of our results to the following NLS type equation with linear dissipation 

idtu + j{t)Au+\u\P^^u + i(Tu = 0, a > 0. 

This equation models wave propagation in dispersion managed fibers including the 
effects of absorption (or damping) by the fiber. In the case without dispersion 
management the influence of the damping term on the possibility of finite time 
blow-up is well studied, cf. [15, 27]. 

Another possibility would be to consider models with only partial dispersion 
management, i.e. where the dispersion management only appears in one space 
direction. An example of this sort can be found in [1] , where the authors study the 
following equation (in d ~ 2): 

idtu + {^{t)dxx + dyy)u + = 0. 

Note that with our choice of 7(t) this equation periodically switches in between the 
usual, "elliptic" NLS and a non-elliptic, or hyperbolic, NLS (in the terminology of 
[18] and [32], respectively). To this end, we recall that the Cauchy problem for the 
non-elliptic NLS 

idtu — dxxU + dyyU + \u\'^u = 

is locally well-posed for initial data in i7^(R^), see [18]. It is conjectured that this 
local solution is actually global. One of the main differences between elliptic and 
non-elliptic NLS is the absence of localized solitary waves for the latter. 

Acknowledgment. The third author wants to thank I. Nenciu for helpful 
discussions. 
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